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Abstract 

We present a simple analysis of the force noise associated with the mechanical damping of the motion of 
a test body surrounded by a large volume of rarefied gas. The calculation is performed considering the 
momentum imparted by inelastic collisions against the sides of a cubic test mass, and for other geometries 
for which the force noise could be an experimental limitation. In addition to arriving at an accurate estimate, 
by two alternative methods, we discuss the limits of the applicability of this analysis to realistic experimental 
configurations in which a test body is surrounded by residual gas inside an enclosure that is only slightly 
larger than the test body itself. 
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1. Introduction 

Gas damping of the motion of a macroscopic test 
body is a potentially sensitivity-limiting source of 
Brownian noise in variety of experiments that are 
sensitive to very small forces. Gas damping in the 
molecular flow regime is characterized by a viscous 
damping coefficient f3 tr = 

firot : 

merous torsion pendulum experiments to be propor 
tional to the residual gas pressure p [H, 0, 0|. The 
power spectrum of Brownian force noise associated 
with the molecular impacts is related to the damp- 
ing coefficient via the fluctuation-dissipation theo- 
rem, which gives Sf — AkBTf3t r - A quantitative 
study of the gas damping and consequent force noise 
acting on small levitated spheres was performed by 



4I7 in translation or 

dV 

in rotation, and has been found in nu- 
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Hinkle and Kendall 0, 0] and experimentally verified 
the Brownian character of residual gas force noise. 

The order of magnitude of the translational mo- 
tion gas damping coefficient can be found easily as 
follows: a test body - or test mass, referred to here 
as TM - with section A that is moving with veloc- 
ity of magnitude V in a gas with molecular number 
density n (related to pressure and temperature by 
n = jrjpf) will be struck, on average, by order nAV 
more molecules per second on the "upwind" side than 
on the "downwind" side. Each molecule imparts a 
momentum of order wiqVt where mo is the mass of 
the molecule and vr = \ / is the characteristic 

1 V mo 

thermal velocity. This gives a translational damping 
coefficient fitr ~ nA x m^VT ~ 

This general dependence is found, albeit with dif- 
ferent prefactors, in the analysis of gas damping force 
noise from damping of mirror vibrational modes for 
gravitational wave interferometers and other oscilla- 
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tors [1, 0, [H , and for the force noise on the cubic TM 
for a space gravitational wave interferometer [§]. To 
get an exact number, we must consider not only the 
momentum exchange normal to the surface, but also 
the role of shear forces acting parallel to the TM sur- 
face. Additionally, we must consider the test mass 
recoil from molecules leaving the surface, which is 
a process that is correlated, via the conservation of 
particles, with the incoming collisions. 

This article presents a calculation of the gas damp- 
ing coefficient, with particular attention to the cubic 
TM geometry. This has been motivated by the need 
for an accurate modeling of gas damping noise on the 
cubic TM that serve as geodetic reference masses in 
the Laser Interferometry Space Antenna (LISA [Io| ) 
and LISA Pathfinder (LPF [11]) missions, in which 
the TM must be free of spurious accelerations at the 
level of fm/s 2 /Hz 1//2 . The dissipation analysis for a 
cube is demonstrated first with a force noise evalua- 
tion and then confirmed with direct calculation of the 
damping coefficient. The force noise analysis is also 
extended to other simple geometries for both trans- 
lation and rotation. 

The calculation presented is straightforward and 
based on the simple physics of diffuse scattering, 
used, for instance, in estimating molecular flow con- 
ductances in tubes. It is of current value, however, 
for several reasons. First, a quantitative calculation 
of gas damping is important for the design of gravita- 
tional wave experiments, and previous estimates ap- 
pear to have underestimated this effect by roughly an 
order of magnitude in noise power [9| . The model pre- 
sented is also easily extended to, and can be tested by, 
torsion pendulum small force experiments. Finally, 
and most importantly, the limitations of the applica- 
bility of this calculation - namely the assumption of 
having a TM surrounded by an infinite gas volume, 
which is discussed in the conclusion - indicate that 
the force noise from gas damping can be grossly un- 
derestimated in many experiments, including LISA, 
which is the subject of a more detailed experimental 
and numerical study [13 ] . The accurate calculation 
of damping in an infinite gas volume is thus the nec- 
essary starting point for studying damping in tighter 
geometries. 




Figure 1: Illustration of inelastic molecular collision with a 
test body surface, with incoming velocity Vi and subsequent 
rcemission with outgoing velocity v . Random, diffuse scat- 
tering, reemission with a cos O distribution is assumed in the 
calculation. 

2. Calculation 

2.1. Force noise calculation 

2.1.1. Force noise normal and parallel to a surface 
element 

We can calculate the force noise spectrum per unit 
area - which we define S± and S\\ for components 
normal and parallel to the surface - from the mean 
square fluctuations of the time average force exerted 
in a time interval Tq. For uncorrelated molecular 
impacts, which we treat as delta- functions in time, 
the time average force on a surface element AA can 
be written 

2 [To ^ 

Fj_ = — dtF± (t) = — y~] m Q {v lk - v ok ) .u± 
lo Jo J o Y 

(1) 

Here, Vik and v k are the incoming and outgoing ve- 
locities of the molecule involved in collision k (an 
analogous expression holds for the force components 
parallel to the surface). 
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For white force noise, as would be expected for mo- 
mentum transfer from uncorrelated, delta-function 
collisions of single molecules against a surface in the 
regime of molecular flow, the single-side noise spec- 
trum can be expressed 



(2) 



S ± AA = 2T ((F ± -(F ± )) 



with an analogous expression for S\i , the force noise 
per unit area parallel to a surface element. 

The calculation can be simplified by considering 
a surface element small enough - and time inter- 
val short enough - such that the probability of more 
than one collision becomes negligible. In this limit - 
uvt TqAA cl - the calculations of the expectation 
value of the average force and its mean square be- 
come simple integrals over the probability, P (vi, v a ), 
of a single collision occuring, in the interval To, with 
incoming and outgoing velocities Hi and v a , 



(F x ) = ^ fd 3 Vid 3 v P(vi,v ) (Vi - v ) .fix (3) 
Jo J 

<-Fx> = d s v i( i 3 v P (vi,v ) [(Vi - v ) .u ± } 2 (4) 
J o J 

The probability P (vi, v a ) can be expressed 

P (vu v ) d 3 Vid 3 v = 

1 \ 3/2 v 2 
n i - — - ) exp--^- (AAcos9iViT ) (vjdvidSli) 



2nv^ 



2vi 



x - — r cos o v exp - — \ (v 2 dv dVt ) 
2irv^ 2v& 



angle independent of the incoming angle, and ree- 
mission is assumed to occur immediately after the 
impact. Weighting by the factor v o cos0 o , which en- 
sures that the incoming and outgoing distributions 
are identical, is the cosine law for diffuse scattering, 
which is used, with substantial experimental verifi- 
cation, for estimating Knudsen diffusion in pipes and 



n pipes anc 
H EI El) 



other geometries (see, for instance, Refs. 
We will discuss the impact of other assumptions of 
the molecule - TM interaction at the conclusion of 
the calculation. 

We can calculate the expected average force, with 
(ili — v ) .u± = Vi cos 9i + v Q cos 9 and restriction of 
the integral over the incoming and outgoing veloci- 
ties to the positive hemisphere o in the range 
[<W2]): 



(Fx) = moAAnvr = pAA 



(6) 



with the final obvious result emerging with the ideal 
gas law, p = nksT. 

In similar fashion we can integrate to obtain the 
mean square average force, 



(7) 



Inspection of Eqns. [6] and [7] shows that \F%j 3> 
(F±) 2 in the studied single-collision regime, and 
Eqn. [2] simplifies to S± = (f± 2 ), and thus 



so 



Si 



(5) 



1/2 



/ 32m Q k B T , 



(8) 



where the velocities are expressed with magnitudes 
v and solid angles Q (angles <f> and 0). The first 
term in Eqn. [5] is the probability of having a 
molecule, assumed to obey a Maxwell-Boltzmann dis- 
tribution with number density n, in the tube of vol- 
ume dA cos 9i Vi T that will hit the wall element in 
the next time interval T with velocity v\. The second 
term corresponds to the probability, given the arrival 
of a molecule, of reemission with speed v into solid 
angle tt . The impacts are assumed to be completely 
inelastic and "memory erasing," with the outgoing 



We note that the "extra" contribution ? comes 
from the cross term (vi cos 0iV o cos0 o ), which can be 
thought of as a correlation between the momentum 
imparted normal to the surface, which is always pos- 
itive, by the arrival and subsequent reemission of a 
molecule. The remaining term - unity in parentheses 
in Eqn. [S] - comes from equal contributions of the 
incoming and outgoing molecules. 

If, in alternative, we were to consider the incom- 
ing and outgoing particle fluxes as uncorrelated pro- 
cesses, which would be true for molecules that stick 
to the surface for a long duration - much longer than 
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the reciprocal of the frequency where we consider the 
force noise - before reemission, the contribution ? 
disappears. Another interesting case is that of elas- 
tic collisions with specular reflection, which is equiv- 
alent to replacing the outgoing distribution with a 
delta function for the specular reflection of the in- 
coming velocity, the factor (l + f J is replaced by a 
factor 2. 

We can perform a similar calculation for either of 
the force components parallel to the surface. For this 
we replace, in evaluating Eqn. 01 the momentum ex- 
change normal to the surface with one of the compo- 
nents parallel to the surface 0, 

wo (vi — v ) = Vi sin 6^ cos <pi — v a sin6* cos <p 

(9) 

For this parallel component, the expectation average 
force is zero. Also, in the mean square of the av- 
erage force, there is no effect of correlation between 
the incoming and outgoing molecules, as the azimuth 
angles are independent of one another and both have 
zero mean momentum exchange parallel to the sur- 
face. We thus obtain 



S 



P 



1/2 



(10) 



Considering the alternative cases discussed above, 
there is no difference in the case of sticking molecules, 
as the incoming and outgoing processes are uncorre- 
cted. In the case of specular reflection, the velocity 
parallel to the surface element is conserved and thus 
S\\ = 0. 

Finally, we note that in the diffuse scattering hy- 
pothesis there is no cross-correlation between the dif- 
ferent components of the force noise, given the az- 
imuthally random reemission of molecules. 

2.1.2. Total force and torque noise on a cubic TM 
and other geometries 
The force noise normal (Eqn. [5]) and parallel (Eqn. 
[TOf to a unit surface element can be integrated over 
the surface of a test body to give the total force along 



1 The resulting noise is the same for the other parallel com- 
ponent, with sin<j> replacing cos</> in the calculation 



a chosen axis. For a cubic test mass, this calculation 
is particularly simple, with the relevant force noise 
component per area (S z for the z axis) uniform on 
each face, given by S± on the two faces (Z) normal 
to z and by Sn on the four faces (X and Y) parallel 
to z. For a cube of side length s, we obtain a total 
force noise 



Sp = 2 dxdyS± + 2 / dydzS\\ + 2 / dxdzS\\ 
J z Jx Jy 

2 f bl2m k B T \ 1/2 f ^ 

We can also evaluate the torque noise along a cho- 
sen axis, which is of relevance to torsion pendulum 
experiments. Considering the torque along the z axis, 
N z = xF y — yF x , we integrate the torque noise per 
area for every surface element, Sn z — x 2 S y + y 2 S x , 
over the faces of the cube. For the Z faces, only shear 
forces, in x and y contribute to the torque. The four 
lateral faces each have contributions from both shear, 
with armlength |, and normal forces, with an arm- 
length that varies across the cubic face. Integrating, 
we obtain 



Sn — 



2 J dydz y 2 S^ + (^) 2 S 
+2 J dxdz x 2 S ± + (0 2 S| 
+2 / dxdy [x 2 5|, +y 2 S||] 



s 4 4s 4 



1/2 



4 / 32 m k B T \ ■ f 7r 



(12) 



The calculation of the total force and torque noise 
has also been performed for a sphere and a right 
cylinder, by straightforward integration of Eqns. [S] 
and [TU] over the TM surface (results in Table 1). We 
note that the rotational damping (or, equivalently, 
the torque noise) on a sphere or cylinder, the latter 
along its axis of symmetry, is due entirely to the shear 
forces and thus would be zero in a model with elastic 
molecular collisions. 
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shape 



Ptr/ [P/V T ] Prot/ \p/v T ] 



f) 1/2 (l+f) 



,1/2 
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^ 2 (if) 1/2 (l + f) 



^ 4 (I) 1/2 



o 



^ 2 (f) 1/2 (l+^+f) 



^ 4 (^) 1/2 (l+f 




7T R 



^) 1/2 x 



i? 
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Table 1: Calculated translational and rotational damping coefficients for different test bodies with cubic, spherical, and cylin- 
drical shapes, along the axes indicated in the figures at left. Note that the damping coefficients are related to the force and 
torque noise spectra via the fluctuation-dissipation formulas, Sp = 4kgT0tr and Sjv = ^kgT0 ro t- 



We note that the results of Table 1 have been con- 
firmed to roughly the percent level, for these three 
geometries and in both translation and rotation, by 
a numerical simulation that traces the impulses of a 
single molecule as it impacts the test mass inside a 
much larger volume, with the reemission distribution 
governed by the diffuse scattering statistics of Eqn. 
[5] [lq. The simulation technique is described in Ref. 
121. 



2.2. Direct calculation of viscous damping coefficient 
for a cube 

As a simple crosscheck for the force noise calcu- 
lated for a cube in the previous section, we can di- 
rectly calculate the viscous damping coefficient f3 by 
calculating the force on a cubic test body moving 
with velocity V — V±u± with respect to the sur- 
rounding gas. Damping will occur both due to the 
molecules which hit the test mass on all faces - with 
a velocity distribution that becomes asymmetric, in 
the test mass reference frame, due to V± - and due 
to the reemitted molecules - which, though emitted 
with the same velocity distribution on all sides of the 
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test mass, still produce a net force, as more molecules 
will be emitted from the upwind side. 

The expectation average force on a surface element 
on the upwind or downwind faces - the faces perpen- 
dicular to the test mass motion - can be calculated, 
as in Eqn. [3J by integrating the momentum exchange 
mo {vi± — Vo±) over the joint probability function, 
Pv ± , which assumes an asymmetry in the reference 
frame of the moving test mass, 

P v± {vi,v ) d 3 Vid 3 v = 

z ' 2 (vix + Vxf + v^ 
exp 



2«| 



x (AAv l± T Q ) (vidvidQi) 
v 2 

v oJ _ exp -— \ {v 2 dv d$l ) 



(13) 



The integrals can be simplified in the limit of 
V± <C vt and neglecting terms of order V 2 , to yield, 
respectively, for the downwind and upwind faces, 



±1 - Vj 



2m 



1/2 



1 



4/ 



(14) 

Again, had we assumed elastic collisions, which essen- 
tially multiplies the force of the incoming molecules 
by two, the factor (l + would become 2, which 
would produce the same result obtained by Christian 
8] for the difference in pressure on two sides of a 
membrane. 

We can similarly calculate the expected force com- 
ponent for a surface element on the faces parallel to 
the TM motion, following Eqn. [9j which yields 



(F\\) = -pAAV\\ 



m 



1/2 



(15) 



Equations [T5] and [13] can be integrated - actually 
a simple multiplication by the surface area s 2 - over 
the six faces of the cube. The leading term in Eqn. [T?] 
cancels out upon summing the upwind and downwind 
faces, to yield a pure velocity damping force, 

^-^= ps *(^) 1/2 (l+*) (16) 
H dV x F \irk B T ) V 8/ K ' 



This agrees with the value of /3, given in Table 1 
and calculated from the force fluctuations and thus, 
as expected, agrees with the fluctuation-dissipation 
theorem. 



3. Discussion 

The expressions derived here for the gas damping 
on a cubic test mass are thus confirmed by two inde- 
pendent calculations. Both calculations of course rest 
on the assumptions of diffuse scattering of molecules 
along the test mass surface. This would be fairly easy 
to test, for the cube (or similar plate) and cylinder 
geometries for the pressure dependence of the rota- 
tional damping coefficient in a torsion pendulum. We 
also note that measuring the rotational damping for 
a cylinder is a direct test of the diffuse scattering hy- 
pothesis, as elastic scattering would give no rotational 
damping. 

Applying Eqn. [IT] to a 46 mm cubic TM, as cur- 
rently employed for LISA, the expected force noise 
from molecular impacts in an infinite gas volume 
would be 



S]l 2 = 0.75fN/Hz 



1/2 



1 ^Pa 



n 



46 mm 



m 
30 m. 



1/4 



(17) 



Two points are important here. The first is the 
magnitude of this term, which, for a given pressure 
and test mass dimension, is roughly three times larger 
than that obtained in a previous analysis Given 



1/2 



1 II 

• p', more am- 



the slow scaling with pressure, u F 
bitious gravitational wave missions, such as BBO 17] 
and DECIGO [l8j], which require lower acceleration 
noise, will have to reach much lower pressures than 
the 10 -6 Pa value quoted here, or employ much larger 
test masses, with a larger mass more than offsetting 
the larger surface area in Eqn. 1171 In general, for 
given TM density, the resulting acceleration noise is 
given by Sl /2 - p 1 / 2 S - 2 - p 1 / 2 M~ 2 ^ for a TM of 
mass M and dimension s. 

The second consideration concerns the applicabil- 
ity of this calculation to frequently encountered ex- 
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pcrimental conditions in which the test mass is sur- 
rounded by an enclosure which is only slightly larger 
than the test mass itself, or more generally is very 
close to the surrounding apparatus. This is certainly 
true for LISA, where the 46 mm cubic test mass is 
enclosed, with a gap of several mm, by the walls of 
a surrounding electrostatic position sensor [l9[ . This 
is also true for terrestrial gravitational wave obser- 
vatories such as Advanced LIGO {2(j, where sus- 
pended TM are in close proximity to similar sus- 
pended masses used for actuation and thermal com- 
pensation. Similar physics of gas damping in a re- 
stricted geometry is fundamental to recent studies of 
"squeeze damping" in MEMS resonators [2 1|. 

For such cases, the assumption, used in the calcu- 
lations presented here, of being able to calculate the 
force noise from the probability of a single impact - 
which is reasonable only if that impact is uncorre- 
lated with previous and future impacts of the same 
molecule - collapses. In a tightly restricted geome- 
try, a molecule emitted from a TM surface will, with 
high probability, strike the opposing enclosure wall 
only to then return to strike the same TM surface, 
imparting momentum with the same sign as the pre- 
vious impact. Random walking, by molecular diffu- 
sion, from one side of the test mass to the other, will 
require many collisions and create a grouping of cor- 
related collisions with force impulses with the same 
sign. This increases force noise, as the momentum 
imparted does not average out as quickly as for un- 
correlated collisions in an infinite gas volume. 

The random walk process down the narrow chan- 
nels between test mass and enclosure is the same that 
determines the flow impedance in a tube. Given that 
test mass motion, inside a tight enclosure, will require 
gas flow around the test mass, this impedance cre- 
ates a pressure head, thus giving an alternate picture 
of the velocity-dependent force, or dissipation, that 
must accompany the increased force noise. These 
constrained volume effects are studied both exper- 
imentally and with numerical simulations, in a re- 
cently published parallel study [12j, which, in addi- 
tion to obtaining the force and torque noise formulas 
derived here in the limit of very large volume, demon- 
strates significantly increased damping in the case of 
gaps that are smaller than the TM itself. 
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